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ABSTRACT 

Context. The origin and heating mechanisms of warm (50< T <200 K) molecular gas in low-mass young stellar objects (YSOs) are 
strongly debated. Both passive heating of the inner collapsing envelope by the protostellar luminosity as well as active heating by 
shocks and by UV associated with the outflows or accretion have been proposed. Most data so far have focussed on the colder gas 
component. 

Aims. We aim to characterize the warm gas within protosteller objects, and disentangle contributions from the (inner) envelope, 
bipolar outflows and the quiescent cloud. 

Methods. High-y CO maps ('"CO 7=6-5 and 7-6) of the immediate surroundings (up to 10,000 AU) of eight low-mass YSOs are 
obtained with the CHAMP+ 650/850 GHz array receiver mounted on the APEX telescope. In addition, isotopologue observations of 
the "CO J =6-5 transition and [C I] Pi line were taken. 

Results. Strong quiescent narrow-line '^CO 6-5 and 7-6 emission is seen toward all protostars. In the case of HH 46 and Ced 110 
IRS 4, the on-source emission originates in material heated by UV photons scattered in the outflow cavity and not just by passive 
heating in the inner envelope. Warm quiescent gas is also present along the outflows, heated by UV photons from shocks. This is 
clearly evident in BHR 71 for which quiescent emission becomes stronger at more distant outflow positions. Shock-heated warm gas 
is only detected for Class flows and the more massive Class I sources such as HH 46. Outflow temperatures, estimated from the CO 
6-5 and 3-2 line wings, are ~100 K, close to model predictions, with the exception of the L 1551 IRS 5 and IRAS 12496-7650, for 
which temperatures <50 K are found. 

Conclusions. APEX-CHAMP+ is uniquely suited to directly probe the protostar's feedback on its accreting envelope gas in terms of 
heating, photodissociation, and outflow dispersal by mapping ~ 1' x 1' regions in high-y CO and [C I] lines. Photon-heating of the 
surrounding gas may prevent further collapse and limit stellar growth. 



X 1- Introduction TMC-l dCabrit & Bertoutlll992l lHogerheijde et al."l998'). Such 

H , molecular outflows are observed to widen with time an d influ - 

P. . Low-mass (M < 3 M^) Stage and Stage 1 protostars ence the evolution of the envelope (e.g. JArce & Sargenjl2005h . 
JRobitaille et alJ200ol) have two distinct components. First, they 

are surrounded by a protostellar envelope, consisting of a reser- Envelope models can be constructed using either a 1-D 

voir of gas and dust, which can feed the star and disk. The enve- or 2-D self-consistent dust radiative transf er calculation con- 

lope can be divided into cold (T ~ 10-20 K) outer regions and a strained by observatio ns of the cold dust (e.g. lShirley et al.l2002l 

warm (T > 50 K) inner region (Adams et al.'l987,'Shirlev et al.' ' J0rgenseneLal.' '2002") and/or the Spectral Energy Distribution 

'2000, Schoier et al. 2002, Shirley et al. 2002, J0r2ensen et al. (S EP) (e.g. W hitney et al. 2003b: a, Robitaille et al. 200l 

20021 IWhitnev et al.1 l2003al iRobitaiUe et al.1 |2006|) . Second, iCrapsi et alJ2008l) . Combining heterodyne observations with ra- 



most, if not all, embedded protostars have molecular out- diative transfer calculations, one can in turn analyze line emis- 
flows (Bachiller 1996, Bachiller & Tafalla 1999, lArce & Sargent! sion of a large range of molecules and determine their origin 



120051, iLee & Hoil2005i, iHirano et al.ll2006) . These outflows have and excitation conditions, both in the inner and outer regions of 

been observed to have a wide variety in intensity, collima- protostellar envelopes as well as in circii mstellar disks and fore- 

tion and aff'ected area, ranging from the strong and large (ten ground components (e.g. lBoogert.MjlJ[2Q 02, J0rg ensenll2OO4 

fliousands of AU) outflow s seen in L 1448, L 1551 IRS 5 iMaret et alJl2 004 lBrinch e t alJ l2007l lL ee et al.,.2007i) . Envelope 

and HH 46 (e.g., Bachille r et al.1 11994 iHeathcote et al.1 1 19961 models predict that flie inner regions are warm (T > 100 K), 

Ivan Kempen et al.|[2009 . from here on referred to as Paper I) to dense (n(H2) > 10^ cm"^) and relatively smaU (R < 500 AU). 

the weak, much more compact, outflows as found for TMR 1 and Most molecular lines observable at (sub)mm wavelengths trace 

the colder gas, but a few lines, such as the high-/ CO lines, di- 

Send offprint requests to: T. A. van Kempen reedy probe the warm gas. The atmospheric windows at 650 
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and 850 GHz are the highest frequency windows in which ob- 
servations of CO (up to energy levels of ~150 K) can be rou- 
tinely carried out. Unfortunately, few studies have succesfully 
observed CO transitions or its isotopologues in these atmo- 
spheric windows due to the excellent weather conditions nec- 
essary. In addition, such studies are o ften limited to single spec- 



Table 1. Adopted Champ^ settings 



tra of only a handful of YSOs (e.g. Schuster et al. 1993; 1995 
HogerMide et al. 1998, CeccarelU et al. 2002, Stark et al. 2004 
Parise et all l2006l Ivan Kempen et al.l l2006l) . The lack of spa 



tial information on the warm gas distribution has prevented 
an in-depth analysis. Comparison between ground-based high- 
J CO observations and far-IR CO transitions (CO 15-14 and 
higher) do not always agr ee on the origin of the high-7 CO lines 
(I van Kempen et al.ll2006h . 

Complex molecules such as H2CO and CH3OH, emit in 
high excitation lines at longer wavelengths, and have been 
observed to hav e surprisingly high abundances in low-mass 
protostars (e.g. Jyan Disho eck et all '1995', 'Blake et al. '1995', 
Ceccarelli etdl 12000, Schoier et alj i200Z Maret et al.. .2004^ 



J0rgensen 2004, Bottin elli et al.ll2004tl2007l) . Unfortunately, the 



abundances of these molecules are influenced by the gas-phase 
and grain surface chemistry, complicating their u se as tracers 
of the physical structure (e.g.. lBisschop et aTlbOOTb . It has been 
proposed that the emission of such molecules originates inside a 
hot core region, a chemically active area close to the star, coin- 
ciding with the passively heated warm inner regio n of the pro- 
tostellar envelope where ices hav e evaporated (C eccarelli et al.l 
l2000llBottinem et al. 2004; 2007h . However, to fully understand 
the origin of these complex organics, knowledge of the structure 
of the warm gas is an essential ingredient. 

Warm gas near protostars can have different origins than pas- 
sive heating alone. Outflow shocks passing through the enve- 
lope can be a source of heat. Quiescent gas, heate d by X-rays or 
UV, is present in the inner envelope. For example. IStauber et al.l 
(12004') show that significant amounts of Far Ultra- Violet (FUV) 
or X-ray photons are necessary to reproduce the line intensi- 
ties and the derived abundances of molecules, such as CO^, CN, 
CH^ a nd NO in both high and low-mass protostars. ISpaans et al.l 
( Il995h investigated photon heating of outflow cavity walls to re- 
produce the observed line intensties and widths of '^CO 6-5 and 
'^CO 6-5 emission in Class I sources (Hogerheiide et al. 1998). 
In paper I, an extension of this model is proposed in which a 
Photon Dominated Region (PDR) at the outflow/envelope cavity 
walls is present in the HH 46 outflow to explain the relatively 
strong and quiescent high-/ CO emission. The emission of [C I] 
can constrain the color and extent of the more energetic photons. 
Apart from the accretion disk, UV photons are also produced by 
the jet shocks in the outflow cavity and the outflow bow shock. 
All these models provide different predictions for the spatial ex- 
tent of the warm gas (~1' for the outflow to <1" for a passively 
heated envelope), as well as the different integrated intensities 
and line profiles of the high-/ CO transitions. 

The most direct tracers of the warm (50 < T < 200 K) 
gas are thus the high-/ CO lines. So far, studies to directly de- 
tect the warm gas components through these lines have rarely 
been able to disentangle the envelope and outflow contributions. 
Far Infrared (IR) transitions of even higher-7 CO, the far-IR CO 
lines, have been observed using the ISO-LWS ins t rument to trace 
the inner regions (e.g., Ceccarelli et al. 1998, Giannini et ajj 
[1999, Nisini et al. 1999, Giannini et al. 2001, Nisini et al. 2002T), 
but could not unambigiously constrain the origin of the warm gas 
emission, due to the limited spatial and spectral resolution. 



Setting 


Lines 


Freq. (GHz) 


Mode 


A 


'-CO 6-5 


691.4730 


OTF/Hexa 


A 


'-CO 7-6 


806.6652 


OTF/Hexa 


B 


'^CO 6-5 


661.1067 


Hexa 


B 


[C I] 2-1 


809.3435 


Hexa 


C 


C'*0 6-5 


658.5533 


Stare 


C 


''CO 8-7 


881.2729 


Stare 



The Chajnantor plateau in northern Chile, where the recently 
commissioned Atacama Pathfinder Experiment (APEXfl is lo- 
cated, currently is the only site able to perform routine obser- 
vations within the high frequency atmospheric windows at high 
(<10") spatial resolutions. The CHAMP^ instrument, developed 
by the MPIfR and SRON Groningen, is the only instrument in 
the world able to simultaneously observe molecular line emis- 
sion in the 650 and 850 GHz atmospheric windows on sub- 
arcminute spatial scales and is thus ideally suited to probe the 
warm gas directly through observations of the 6-5, 7-6 and 8-7 
transitio ns of CO and its isotopologu es with 7-9" angular res- 
olution ( Kasemann etal.1 12006. Gust en et al.ll2008l) . CHAMP+ 
has 14 pixels (7 in each frequency window) and is thus capa- 
ble of fast mapping of the immediate surroundings of embedded 
YSOs. The Herschel Space Observatory will allow observations 
of far-IR CO lines at spectral and spatial resolution similar to 
APEX (~10"). Note also that the beam of Herschel is compara- 
ble or smaller than the field of view of CHAMP+ at its longer 
wavelenghts (~500//m). Thus the CHAMP+ data obtained here 
provide information of the distribution of warm gas within the 
Herschel beams. 

In paper I, we presented the results for one source, HH 
46 IRS. In this paper, we present observations of CO and its 
isotopologues using CHAMP+ for seven additional embedded 
YSOs and compare them with the HH 46 case. §2 presents the 
sample and observations. The resulting spectral maps are shown 
in §3. §4 discusses the envelope and outflow structure, while the 
heating within protostellar envelopes and molecular outflows is 
analyzed in §5. §6 investigates the relation between the emis- 
sion of more complex molecules and the emission of high-7 CO. 
Final conclusions are given in §7. 



2. Sample and observations 

2.1. Observations 

The sample w a s obs er ved using th e CH AMP+ array 
(Kasemann et al.' '20061 iGiisten et all |2008|) on APEX 
(Kasemann et al. 2006). CHAMP"'' observes simultaneously 
in the 650 GHz (450 yum, CHAMP+-I) and 850 GHz (350 
fim, CHAMP^-II) atmospheric windows. The array has 7 
pixels for each frequency, arranged in a hexagon of 6 pixels 
around 1 central pixel, for a total of 14 pixels. During the 
observations, the backend consisted of 2 Fast Fourier Transform 
Spectrometer (EFTS) units serving the central pixel, and 12 
MPI- Auto-Correlator Spectrometer (MACS) units serving the 
other pixels. The EFTS units are capable of observing up to a 
resolution of 0.04 km s"' (0.12 MHz) and the MACS units to a 



' This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fiir Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 
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Source 


RA 

(J2000) 


Dec 
(J2000) 


D 

(pc) 


^^bol 

(io) 


Tbo\ 
(K) 


Class 


NGC 1333 IRAS 2 


03:28:55.2 


+31:14:35 


250 


12.7 


62 





L1551 IRS 5 


04:31:34.1 


+ 18:08:05.0 


160 


20 


75 


I 


TMR 1 


04:39:13.7 


+25:53:21 


140 


3.1 


133 


I 


HH46 


08:25:43.8 


-51:00:35.6 


450 


16 


102 


I 


Ced 110 IRS 4 


11:06:47.0 


-77:22:32.4 


130 


0.8 


55 


I 


BHR71 


12:01:36.3 


-65:08:44 


200 


11 


60 





IRAS 12496-7650 


12:53:17.2 


-77:07:10.6 


250 


24 


326'' 


I 


RCrA IRS VA"^ 


19:01:55.2 


-36:57:21.0 


170 


- 


- 







Sett." 


Ref. Cont. 




Ref.CO 






NGC 1333 IRAS 2 


A 


1,2,3,4,5,6,7 




vii 






L1551 IRS 5 


A 


2,8,9,10,11,12 




i 






TMR 1 


AB 


13,14,15,16 




i 






HH46 


ABC 


1,17,18 




ii,iii 






Ced 110 IRS 4 


AB 


7,19,20,21,22 




iii,iv 






BHR71 


AB 


17,23,24,25 




V 






IRAS 12496-7650 


AB 


26, this work 




iii,viii 






RCrA IRS 7A'^ 


A 


1,7,27,28 




iii, viii 







" Observed settings with CHAMP+. See Table[T] 

* IRAS 12496-7650 is likely to be viewed face-on and thus has strong IR emission and a high Tboi- See van Kempen et al. 2009 submitted 

RCrA IRS 7 is a binary. Pointin g was chosen to be on RCrA IRS 7A, believed to be the main embedded source, following results from van 
Kempen et al. 2009 submitted and lGroppi et al.l ( 120071) . RCrA IRS 7B is located 15" to the east. At mid- and far-IR wavelengths, the source is 

heavily confused with RCrA. Therefore, no reliable Tboi and Lboi could be derived. 

Contin uum References: I: Di Francesco et al. (2008) 2: Motte & Andre (2001) 3: lGutermuth et al.l j2008l) 4 : lEnoch et al.l j2006h 5:|Hatchell et al. 



2007h 6:|Sandell & Knee (2001) 7: Froebrich (2005) 8: Chandler & Richer (2000) 9: Reipurth et al. (2002) 10: Osorio et al. (2003) IljLiuetal 
19961) 12: Butneretal. (1991) 13: Kenvon et al. (1994) 14: Hogerheijde et al. (1999) 15: Furlan et al. (2008) 16: Terebey et alJ (|I993|) 17: 
Hennine & Launhardt (1998) 18: van Kempen et al. (2009) 19: Reipurth et al. (1993) 20: Luhman et al. (2008) 21: Lehtinen et al. (2001) 22: 
Lehtinen et al. ( 2003) 23: .Evans et al, t2007.) 24: .Bourke et al.. ( .1997.) 25: Bourka ( .2001.) 26: .Henning et al., ( ^1993.) 27: ,Nutter et ah t2005» 28: 

Groppietal. (2007) 

CO references : i:iHogerh eiide et al.' [T99§) ii: van K empen et al. 2009 sub mitted iii: Paper I iv: lHiramatsu et afl ( l2007f) v: iParise et al.l 1200^ vi: 
lKnee&Sand"eil ( l2000h vii: Giannini et al. (200l]) viii: lGiannini et al" 

I CHI) 



spectral resolution of 0.37 km s"' (1 MHz), both at a frequency 
of 806 GHz. The observations were done during three observing 
runs in June 2007, October-November 2007 and July 2008 in 
three different line settings; see Table [T] Note that in July 2008, 
all pixels were attached to FFTS units and no MACS backends 
were used. For mapping purposes, the array was moved in a 
small hexagonal pattern to provide a fully Nyquist sampled map 
or in an on-the-fly (OTF) mode. The 'hexa' pattern covers a 
region of about 30" x 30". For some sources, a slightly larger 
area was mapped in a small OTF (40" x 40") map in '^CO 
6-5/7-6. Both kinds of maps were regridded to a regular grid 
with standard rebinning algorithms included in the CLASS 
packag^. A binning method that uses equal weight binning was 
adopted to produce contour maps. A pixelsize of l/3rd of the 
beamsize was adopted to create accurate contour maps. CLASS 
was also used as the main reduction package for individual 
lines. At the edge of the maps, noise levels are often higher due 
to the shape of the CHAMP+ array. NGC 1333 IRAS 2, TMR 
1 and RCrA IRS 7 were observed using the hexa mode and 
have smaller covered areas. Due to the different beams at 690 
and 800 GHz, small differences exist between the areas covered 
in CO 6-5 and 7-6. For the C'^^O 6-5 /"CO 8-7 (setting C; 
only done for HH 46), a stare mode was used to increase the 
S/N within the central pixel. A position switch of 900 " was 
used for all settings, except for the stare setting C, which used a 
beam-switching of 90". 



^ CLASS is part of the GILDAS reduction package. See 
|http://www.rram.fr/IRAMFR/GrLDAS for more information. 



Beam efficiencies, derived using observations on planets dur- 
ing each run, are 0.56 for CHAMP+-I and 0.43 for CHAMP+-II. 
These efficiencies were found to vary by less than 10 % over a 
single observing run and between runs. Between pixels the vari- 
ation is similar, all within 10%. Depending on observing mode, a 
given sky position is covered by many (all) pixels, so the average 
(relative) calibration eiTor will be a few per cent. Typical single 
side-band system temperatures are 700 K for CHAMP^-I and 
2100 K for CHAMP+-II. The 12-m APEX dish produces a beam 
of 9" at 650 GHz and 7" at 850 GHz. Pointing was checked 
on various planets and sources and was found to be within 3". 
Calibration of the sources was done using similar observations, 
as well as hot and cold loads. The typical sideband rejection at 
both frequencies was measured to be less than lOdB. At APEX 
10 dB is used as the input to the calibrator (both continuum and 
line. From all these effects, the total calibration error is estimated 
to be ~30%, including the atmospheric model. Note that all tech- 
nical aspects of CHAMP+ will be discussed more extens ively in 
Glisten et al. (to be submitted to A&A, see also Giiste net al.l 

dloosh ). 

The difference in_C07-6 emission of 60% for IRAS 12496- 
7650 in van Kemp en et al.l (2006) and this paper can be ac- 
counted for by the respective calibration and pointing eiTors of 
FLASH and CHAMP+. A fac tor of 2 difference is foun d with 
the observations of TMR 1 in iHogerheiide et al.l (Il998l) . Since 
a similar factor of 2 is seen between the observations of T Tau 
in lSch uster et al. ( 1995) and Hogerheijde et al. ( 1998), it is most 
likely that the side - band gain ratios for the observations taken by 
IHogerheiide et all (119981) were incorrectly calibrated. 
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Table 3. Properties of the '^CO lines". 



Source 


Offpos. 


Vlsr 




CO 6-5 












/ TundV 




red 


blue 


RMS* 




(",") 


km s"' 


Kkms-' 


K 


Kkms-' 


K 


NGC 1333 IRAS 2 


0,0 


7.0 


34.3 


6.8 


8.9 


13.2 


0.6 




0,20 


7.0 


43.7 


9.8 


20.0 




0.6 




-20,-20 


7.0 


36.3 


7.6 


2.2 


Il.O 


0.7 


L 1551 IRS 5 


0,0 


6.8 


39.5 


7.8 


13.8 


10.7 


0.3 


TMR 1 


0,0 


4.5 


23.4 


8.1 


2.5 


5.2 


0.3 


HH 46 


0,0 


4.4 


45.3 


9.3 


10.6 


10.4 


0.3 




-20,-20 


4.4 


33.4 


9.1 


11. 1 




0.5 


Ced 110IRS4 


0,0 


3.5 


20.8 


7.8 




4.8 


0.4 


BHR71 


0,0 


-4.5 


38.6 


5.2 


9.1 


12.7 


0.4 




10, -40 


-4.5 


96.3 


16.5 




68.1 


I.O 




0, 40 


-4.5 


75.8 


13.0 


41.5 




I.O 


IRAS 12496-7650 


0,0 


1.8 


22.1 


5.9 






0.5 


RCrA IRS VA"^ 


0,0 


5.7 


337.9 


36.3 


33.6 


10.9 


0.3 


Source 


Offpos. 


Vlsr 




CO 7-6 












/ TundV 




red 


blue 


RMS* 




(",") 


km s"' 


Kkms-' 


K 


Kkms-' 


K 


NGC 1333 IRAS 2 


0,0 


7.0 


29.2 


6.3 


8.8 


5.2 


1.7 




0,20 


7.0 


65.8 


lO.I 


27.8 


5.2 


1.7 




-20,-20 


7.0 


- 








3.5 


L 1551 IRS 5 


0,0 


6.5 


52.9 


9.2 


18.9 


11.6 


0.4 


TMR 1 


0,0 


4.5 


19.8 


lO.I 


1.5 


1.8 


0.7 


TUTU /I 

Hrl 46 


0,0 


A A 

4.4 


1 A ^ 


8.4 




9.9 


U.5 




-20,-20 


4.4 


29.8 


6.2 


15.7 




I.O 


Ced 110IRS4 


0,0 


3.5 


23.1 


7.1 






I.O 


BHR71 


0,0 


-4.5 


36.6 


4.4 


5.4 




0.7 




10, -40 


-4.5 


41.5 


7.0 


5.6 


17.2 


2.8 




0, 40 


-4.5 


54.7 


9.1 


27.0 




3.0 


IRAS 12496-7650 


0,0 


1.8 


19.9 


6.1 






1.4 


RCrA IRS 7A^ 


0,0 


5.7 


407.6 


46.3 


29.0 


22.2 


0.7 



" Red is the redshifted outflow lobe calculated from -10 to -1.5 km s"' with respect to the source velocity. Blue is the blueshifted outflow lobe 
calculated from -1-1.5 to +10 km s"' with respect to the source velocity. 
* 1 IT in a 0.8 km s"' channel. 

^ Outflowing gas refers to emission in -20 to -8 km s"' and +8 to +20 km s"' from line center. 



Table 4. CO isotopologue and [C I] properties at the source position. 



Source 


1 


'CO 6-5 




'^CO 8-7 


C'^O 6-5 


[C I] 2-1 




/ TuBdV 


-^peak 


FWHM 


/ TusdV 


jTuBdV 


/ TusdV 




Kkms' 


K 


km s"' 


K km S-' 


K km S-' 


Kkms-' 


TMR 1 


4.3 


1.9 


2.0 






2.9 


HH46 


3.9 


3.1 


1.5 


<0.3'' 


<0.15° 


2.3 


Ced 1I0IRS4 


5.2 


2.1 


2.3 






2.5 


BHR7I 


1.4 


1.5 


0.8 






<0.6" 


IRAS 12496-7650 


3.8 


1.7 


2.1 






<0.6" 



Icr level in a 0.7 km s bin. 



2.2. Sample 

The sample consists of eight well-known and well-studied em- 
bedded protostars, with a slight bias toward the southern sky. 
A variety of protostars in mass, luminosity, evolutionary stage 
and parental cloud is included. All sources have been studied 
in previous surveys of embedded YSOs in the sub-mm (e.g., 
[j0rgensen et al. 2002; 2004, Groppi et al. 20 02). More infor- 
mation on southern sources can be found in Ivan Kempen et alj 
(120061) (IRAS 12496-7650), van Kempen et al. 2009 submitted 



(IRAS 12496-7650, RCrA IRS 7 and Ced 110 IRS 4) and Paper 
I (HH 46). Table |2] gives the parameters of each source and its 
properties. References to previous continuum studies (column 7) 
include infrared (IR) and (sub)millimeter dust continuum pho- 
tometry, and studies of CO emission (column 8) include both 
submillimeter and far-IR lines. Note that the binary source of 
N1333 IRAS 2A, N1333 IRAS 2B is not covered by our map. 
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2.3. Spectral energy distribution 

For all sources, SED information was acquired from the 
literature, ranging from near-IR to the (sub-)mm wave- 
lengths. Spitzer-IRAC (3.6, 4.5, 5.6, 8.0 fj.m), Spitzer- 
MIPS (24, 70 a nd 160 Atm), submillime ter photometry from 
JCMT/SCUBA (I Pi Francesco"etan llOOSi) . SEST/SIMBA and 
IRAM-30m/MAMBO are all included, if available. For sources 
with no reliable Spitzer fluxes, IRAS fluxes at 12, 25, 60 and 
100 jum were used. In addition to the obtained fluxes, the bolo- 
metric luminosity and temper ature, Lboi and Thoi we re calculated 
using the mid-point method jPunham et alj|2008l) . The results 
can be found in Table |2l Note that due to a lack of mid- and 
far-IR fluxes with sufficient spatial resolution, the derived values 
for IRAS 12496-7650 and Ced 1 10 IRS 4 are highly uncertain, 
with expected eiTors of 50%, while no Tboi and Lboi could be de- 
termined for RCrA IRS 7 due to confusion. For other sources, 
errors are on the order of 10-20%. 



3. Results 

3.1. Maps 

Figures[T]to[8]show the integrated intensity and the red- and blue- 
shifted outflow emission mapped in the CO 6-5 line, as well as 
the spectra at the central position of all observed emission lines. 
Both the integrated intensity and outflow maps were spatially 
rebinned to a resolution of 10". Figures |9] to [T2l show the CO 6- 
5 and 7-6 spectra taken in the inner 40" x 40" of all the sources, 
also binned to square 10" x 10" bins. 

Table [3] presents the intensities of both CO lines. Here, the 
total integrated emission, the peak main beam temperature, the 
emission in the red and blue outflow wings shown in the outflow 
maps and the noise levels are given for the central position, as 
well as at positions with clear detections of outflow wings away 
from the center. The limits for the red and blue outflow were cho- 
sen to be - 10 to - 1 .5 km s"' with respect to the source velocity for 
the blueshifted emission and -1-1.5 to +10 km s"' with respect to 
the source velocity for the redshifted emission for most sources. 
These limits were chosen after examination of the profiles of the 
spectra and subtracting gaussians that were fitted to the central 
3 km s ' of the profile. For all sources, except RCrA IRS 7, the 
difference in the blue- and red-shifted emission between the two 
methods was less than 5 %. The outflowing gas of RCrA IRS 7 
was derived by limiting the red- and blue-shifted emission to -20 
to -8 and -1-8 to -1-20 km s"' from line center. This corresponds to 
a FWHM«i8 km s"' for the quiescent material. Broad line widths 
of 3 km s ' are also seen for the rarer isotopologues C"^0 and 
C'^0, much wide r than typica lly observed for these species in 
low-mass YSOs ( Schoier et al.ll2006) . Noise levels (see Table[3]l 
differ greatly between sources and even within single maps. 

'^CO 6-5 and 7-6 was detected at the central position of 
all sources, ranging from 20.8 K km s"' (7'peak=7.8 K) for Ced 
110 IRS 4 for CO 6-5 to 407.6 K km s ' (rpeak=46.3 K) for 
RCrA IRS 7A for the CO 7-6 Une. All maps show extended 
emission, except for IRAS 12496-7650, which shows unresolved 
emission in the CO 7-6 transition. However, the scales on which 
extended emission is seen varies significantly, with detections at 
all mapped positions for RCrA IRS 7 to only 1 or 2 for IRAS 
12496-7650, TMR 1 and Ced 110 IRS 4. AH sources except 
RCrA IRS 7 and NGC 1333 IRAS 2 show spectra with a sin- 
gle peak over the entire map, while the latter two have spectra 
that are self-absorbed. 



Figures l9l to [T2l clearlv identify the variation of the fine pro- 
files across the maps, especially when outflowing gas is present, 
such as in the maps of NGC 1333 IRAS 2, BHR 71, HH 46 and 
RCrA IRS 7. The sources for which little to no shocked emis- 
sion is seen do show spatially resolved CO 6-5 and 7-6 emis- 
sion, but always quiescent narrow emission located close to the 
central pixel. 

3.2. Outflow emission 

From the outflow maps in Fig. [T] to [8] it can be concluded 
that the contributions from shocks within the bipolar outflows 
to the warm gas differ greatly from source to source. RCrA 
IRS 7, NGC 1333 IRAS 2, HH 46, and BHR 71 produce 
spatiafly resolved flows, but TMR 1, L 1551 IRS 5, IRAS 
12496-7650 and Ced 110 IRS 4 do not have warm shocked 
gas that results in broad high-/ CO in any of the off-positions. 
Small shocks on the source position are seen, but are gen- 
erally weak. For all sources in the sample, outflow emission 
has been detected in low (7up < 3) excitation CO lines, 
although such flows have large differences in spatial scales, 
ranging from tens of arcminutes to the central twenty arc 



seconds (e.g. Moriartv-Schieven & Snelll Il988l Ifiid lund et al 
1989, Cabrit & Bertout 1992, Bachill er et al.ll994 lBourke et al 
il997UHog erheiide et al. 1998, Pa rise et al.l 2006', van Kempen et 
al. 2009 submitted). Table |3] gives the integrated intensities of 
the CO emission at selected off-positions for the outflowing gas. 

3.3. Isotopologue observations at the central position 

For HH 46, transitions of ^^CO 7=6-5 and 8-7 and C"*0 7=6-5 
as well as [C I] 2-1 were observed (see Table [T] and |2]l. TMR 1, 
Ced 110 IRS 4, IRAS 12496-7650 and BHR 71 were observed 
in '^CO 6-5 and [C I] 2-1. The results at the central position can 
be found in Fig. [T] to [8] as weU as Table g] All '^CO 6-5 spectra 
can be fitted with single gaussians. However, the width of the 
gaussians varies with a FWHM of ~1.2 km s ' for HH 46 and 
BHR 71 to -2.2 km s ' for Ced 110 IRS 4 and IRAS 12496- 
7650. [C I] 2-1 is detected for TMR 1, HH 46 and Ced 1 10 IRS 
4 (~3 cr). No line was found for IRAS 12496-7650 and BHR 71 
down to a Icr limit of 0.6 K in 0.7 km s ' bins. Integrated line 
strengths are on the order of 3 K km s ' with widths of 0.75 km 
s ' . '^CO 8-7 and C"*0 6-5 were observed towards HH 46 only, 
but not detected down to a Icr level of 0.15 and 0.3 in a 0.7 km 
s ' channel. 



3.4. ^^CO 6-5 and [C I] 2-1 maps 

For '^CO 6-5 and [C I] 2-1 maps were also obtained. Spectra at 
the central position are given in Figs.[T]to|8] Integrated intensity 
maps are presented in Fig. [13] ('^CO 6-5) and Fig. [14] ([C I] 
2-1) It is seen that such observations are often dominated by 
centrally located unresolved emission, but not always peaked at 
the source. Both HH 46 and BHR 71 also show some isotopic 
emission associated with the outflow. For HH 46, '^CO 6-5 is 
only detected ofF-source for the blue outflow, where part of the 
outflow is unobscured by cloud or envelope (Paper I). 



4. Envelope 

4. 1 . Envelope models 

In order to investigate whether high-7 CO emission can be 
reproduced by a passively heated envelope, the properties of 
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Fig. 1. Maps of CO 6-5 (/op row) and CO 7-6 (bottom row) of NGC 1333 IRAS 2. The left-most image shows the total integrated intensity over the line. The middle figures show the 
outflow contributions from the red (dashed lines) and blue (solid lines) outflow. Contours are in increasing levels of 4 K km s~' for both transitions. The outflow contributions are calculated 
by only including emission greater or smaller than +/- 1.5 km s~^ from the central velocity. The right-most image at the top row shows the '"CO spectra at the central position. 
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Fig. 2. Maps of CO 6-5 (top row) and CO 7-6 (hot lorn row) of L 1551 IRS 5 . The left-most image shows the total integrated intensity over the line. The middle figures show the outflow 
contributions from the red (dashed lines) and blue (solid lines) outflow. Contours are in increasing levels of 0.5 K km s"' for both transitions. The outflow contributions are calculated by 
only including emission greater or smaller than +/- 1.5 km s~' from the central velocity. The light-most image at the top row shows the '"CO spectra at the central position. 



the protostellar envelopes were calculated by modelling 850 radiati v e transfer code DUSTY described i n Ivezic & Elitzurl 
fim continuum images, SED information and the 1-D dust (Il997h dSchoier et alJ 120021 |j0rgensen et al] [2002) . Similar to 
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Fig. 3. Maps of CO 6-5 {top row) and CO 7-6 {bottom ravi) of TMR 1 . The left-most image shows the total integrated intensity over the line. The middle figures show the outflow 
contributions from the red {dashed lines) and blue {solid lines) outflow. Contours are in increasing levels of 1 K km s"' for both transitions. The outflow contributions are calculated by only 
including emission greater or smaller than +/- 1.5 km s"' from the central velocity. The right-most image at the top row shows the ^^CO spectra at the central position, while the right-most 
image at the bottom row shows the spectra of the observed isotopologues and [C 1] 2-1 at the central position. 
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Fig. 4. CO 6-5 {left) CO 7-6{middle) maps and spectra of the central position of HH 46. The left-most image shows the total integrated intensity over the line. The middle figures show 
the outflow contributions from the red (dashed lines) and blue (solid lines) outflow. Contours are in increasing levels of 2 K km s"' for CO 6-5 and 3 K km s"' for CO 7-6. The outflow 
contributions are calculated by only including emission greater or smaller than +/- 1.5 km s"' from the central velocity. The right-most image at the top row shows the '-CO spectra at the 
central position, while the right-most image at the bottom row shows the spectra of the observed isotopologues and [C I] 2-1 at the central position (see also paper I). 
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Fig. 5. Maps of CO 6-5 {top row) and CO 7-6 (bottom row) of Ced 110 IRS 4. The left-most image shows the total integrated intensity over the line. The middle figures show the outflow 
contributions from the red (dashed lines) and blue (solid lines) outflow. Contours are in increasing levels of 3 K km s"' for both transitions. The outflow contributions are calculated by only 
including emission greater or smaller than +/- 1.5 km s"' from the central velocity. The right-most image at the top row shows the ^^CO spectra at the central position, while the right-most 
image at the bottom row shows the spectra of the observed isotopologues and [C I] 2-1 at the central position. 
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Fig. 6. Maps of CO 6-5 (top row) and CO 7-6 (bottom row) of BHR 71. The left-most image shows the total integrated intensity over the line. The middle figures show the outflow 
contributions from the red (dashed lines) and blue (solid lines) outflow. Contours are in increasing levels of 10 K km s~' for CO 6-5 and 5 K km s~' for CO 7-6. The outflow contributions 
are calculated by only including emission greater or smaller than +/- 1.5 km s~' from the central velocity. The right-most image at the top row shows the '"CO spectra at the central position, 
while the right-most image at the bottom row shows the spectra of the observed isotopologues and [C I] 2-1 at the central position. 




Fig. 7. Maps of CO 6-5 {lop row) and CO 7-6 {bottom row) of IRAS 12496-7650. The left-most image shows the total integrated intensity over the line. The middle figures show the 
outflow contributions from the red {dashed lines) and blue {solid lines) outflow. Contours are in increasing levels of 3 K km s"' for both transitions. The outflow contributions are calculated 
by only including emission greater or smaller than +/- 1.5 km s"^ from the central velocity. The right-most image at the top row shows the ^^CO spectra at the central position, while the 
right-most image at the bottom row shows the spectra of the observed isotopologues and [C 1] 2-1 at the central position. 



previous studies, the so-calle d OH5 dust opacities were used 
( lOssenkopf & HenningI [19941) . Although NGC 1333 IRAS 2, 
LI 551 IRS 5 and T MR 1 were included in a similar study 
by j 0rgensen et aj] (l2002|). new flu xes from Spitzer (e.g., 
iLuhman et al.l l2008l iGutermuth eTall l2008i see Table ^ have 
since come available and must be included. To constrain the 
models both the SED from ~50 fim to 1.3 mm and the spatial 
distribution of the sub-mm continuum emission at 850 fim were 
used, with the exception of BHR 71 and Ced 110 IRS 4, for 
which no 850 fim map is available. A normalized radial emis- 
sion profile with a power law index of 1 .5 for the 850 fj.m emis- 
sion is assumed for these two sources. The radial profiles of L 
1551 IRS 5, RCrA IRS 7A, NGC 1333 IRAS 2 and TMR 1 were 
obtained from 850 //m irn ages of the processed SCUBA archive 
dPi Francesco etalj|2008h . For HH 46 and IRAS 12496-7650, 
recent data from LABOCA at 870 //m were used (see Paper I 
and Nefs et al. in prep.). 

The parameters of the best-fitting envelope models can be 
found in Table |5] together with the corresponding physical pa- 
rameters of the envelope. The three main parameters of DUSTY, 
Y, the ratio over the inner to outer radius, p, the power law expo- 
nent of the density profile n oc r and Tioo^jm, the opacity at a 
100 fim, are scaled by the Lboi and the distance, D given in Table 
m to get the physical properties of each source. Fig. [TSl shows 
the best-fitting model of each source of the radial profiles of the 
850 jum images and the SEDs of the entire sample. 

The inner radii of the protostellar envelopes range from 5 to 
35 AU corresponding to = 250 K, a limit chosen by us. Most 



sources also show a steep profile with all sources having p > 1 .7, 
with the exception of Ced 110 IRS 4, which has a p-lA. 



4.2. CO emission witliin protostellar envelopes 

Using the best-fit envelope temperature and density structure 
derived from the dust emission, the CO intensities and line 
profiles from the protostellar envelopes were in turn simulated 
with the self-consistent ID molecular line radiative transfer code 
RATRAN (Hogerheiide & van der Tak 2000) using data files 
of the LAMDA database ( Schoier'etaLll2005l). CO abundances 



^ Note that this power law index is not the same as the power law 
index p for the normalized 850 pm radial emission profile 



with respect to H2 are taken from iJOrgensen et"an (l2005h . A 
'drop' abundance with Xo=2.7xlO and X^^IQ'^ is adopted. 
This 'drop' abundance profile describes a warm (T > Tfieeze) 
inner region with a high abundance Xq and a region in which 
T < Tf^eeze and n > «(ie where CO is frozen out to a low abun- 
dance Xj. In the outer region (« < «de), the abundance is again 
high at Xq because freeze-out timescales become longer than 
the typical life-times of protostars. In our models, rfieeze=30 
K and Wde^lO^ cm~^ a re adopted, following the conclusions of 
|j0rgensen et al.l (l2005h . There, the derived abundances are based 
upon the emission of low-excitation optically thin lines such as 
C'**0 lines (both 2-1 and 3-2, 1-0 is often dominated by the very 
cold cloud material) and C'^O. Paper I showed that contrary to 
the low-y CO lines there is little difference between 'jump' and 
'drop' abundances for the emission in high-/ transitions. A static 
velocity field is assumed with a turbulent width of 1 km s ' . Due 
to the static nature of the velocity field, excessive self-absorption 
is seen in line profiles of CO lines up to the 8-7 transition. The 
total area of a gaussian fitted to the line wings is used in those 




Fig. 8. Maps of CO 6-5 (top row) and CO 7-6 (boltom row) of RCrA IRS 7a. The left-most image shows the total integrated intensity over the line. The middle figures show the outflow 
contributions from the red (dashed lines) and blue {solid lines) outflow. Contours are in increasing levels of 10 K km s"' for both transitions. The outflow contributions are calculated by 
only including emission greater or smaller than +/- 8 km from the central velocity. The right-most image at the top row shows the '"CO spectra at the central position. 



Table 5. Inferred envelope properties from the DUSTY models. 



Source 


Y 


P 


TlOO 


^inner 




"lOOOAU 


Menv(< i?outcr) 










(AU) 


(10* AU) 


(lO-* cm-3) 


(Mo) 


RCrA IRS 7A 


700 


2.0 


0.2 


24.7 


5.9 


13 


0.9 


NGC 1333 IRAS 2 


1500 


1.9 


1.8 


31.3 


1.6 


140 


2.2 


L 1551 IRS 5 


500 


1.7 


0.3 


26.5 


1.3 


33 


0.6 


TMR 1 


2400 


1.9 


0.2 


9.5 


2.3 


5.4 


0.1 


IRAS 12496-7650 


2500 


2.2 


0.2 


35.7 


8.9 


8.7 


0.3 


HH46 


700 


1.8 


2.6 


30.0 


2.1 


250 


5.1 


BHR71 


1300 


1.9 


0.2 


15.6 


2.0 


8.5 


0.2 


Ced 110 IRS 4 


1400 


1.4 


0.5 


5 


0.7 


51 


0.6 



cases to derive an upper limit. This is a very strict limit as the 
true CO emission associated with the envelope is best fitted by 
a infall velocity (Schoier et al. 2002), producing integrated in- 
tensities between the two limits. However. ISchoier et al.l (l2002l) 
show that the envelope emission modelled with an infall veloc- 
ity is in the worst case only a factor 2 greater than the intensity 
derived from the static envelope. The lower limit is derived from 
the actual modelled integrated intensity with the static velocity 
field. 

Figure[T6]shows the resulting integrated intensities produced 
by the model protostellar envelopes of all CO lines from 7=1-0 
up to 7=19-18. The data for all transitions were convolved with 
a 10" beam used for CHAMP^, except the three lowest transi- 
tions, which are convolved with a beam of 20". Such beams are 
typical for single-dish submillimeter telescopes. 10" will also 
be the approximate beam for several transitions covered by the 
PACS and HIFI (Band 6 and 7) instruments on Herschel at the 
higher frequencies. Overplotted are the observed line strengths 



from various CO line s of different studies , inc luding the far- 
IR high-y CO lines of ' Giannini et all d 19991) and lGiannini et al.1 
(2001), assuming that all flux observed by the ISO-LWS in its 
80" beam originates in a 10" region. See Table |2] for the refer- 
ences used. 

The low-7 CO emission can often be completely reproduced 
by the envelope models, as it is dominated by the colder gas in 
the outer regions of the envelopes. For HH 46, Ced 110 IRS 4 
and RCrA IRS 7, the observed quiescent emission in the CO 6- 
5 and 7-6 lines is clearly brighter by a factor of 3-5 than the 
modelled envelope emission, even for the '^CO 6-5 lines (see 
Table©. 

For RCrA IRS 7, the difference is almost an order magnitude. 
Deep C'**0 6-5 spectra are needed to fully pin down the enve- 
lope models. In Paper I, the quiescent CO 6 -5 and 7-6 emissio n 
of HH 46 was attributed to 'photon heating' dSpaans et al.ll995h . 
both by UV from the accretion disk and shocks inside the out- 
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Fig.9. Specti-a of '^CO 6-5 over an area of 40" x 40" of NGC 1333 IRAS 2 (/e/r) and L 1551 IRS 5 ((right). NGC 1333 IRAS 2 
was observed using the hexa mode, so the covered area is slightly smaller than that of L 1551 IRS 5 taken in OTF mode. For all 
spectra in both sources, horizontal axes range from -10 to 25 km s"' and vertical axes from -5 to 12 K. For both sources, the main 
outflow axis of the red and blue are shown with a dashed and solid lines. 
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Fig, 10. Spectra of ^^CO 6-5 over an area of 40" x 40" of TMR 1 (left) and HH 46 ((right). TMR 1 was observed using the hexa 
mode, so the covered area is slightly smaller than that of HH 46 taken in OTF mode. For all spectra in both sources, the horizontal 
axes range from -10 to 20 km s"' with vertical axes from -5 to 12 K. For HH 46, the red outflow axis is shown with a dashed line. 



flow cavity. Fig. [16] clearly shows that this method of heating 
likely applies to other sources. 

In contrast, the envelopes of both IRAS 12496-7650 and 
L 1551 IRS 5 can account for all the emission detected 
in the '^CO 6-5 an d 7-6 Hnes. iGiannini etall (fl999 ) and 
iGiannini etal] (1200 ll) report emission of high-7 CO (14-13 to 
19-18) lines at far-IR of NGC 1333 IRAS 2, IRAS 12496-7650 
and RCrA. For RCrA, RCrA IRS 7 is within the beam, but the 
emission is probably dominated by emission from RCrA itself. 
Fluxes in excess of 10"^° W cm"^ are seen for IRAS 12496- 
7650. They assumed that the flux originates within the central 



3" (~400 AU). It is very clear from Fig.[T6]that such observed 
emission cannot be produced by a passively heated envelope. A 
more likely explanation is that the CO emission detected with 
ISO-LWS is either located outsi de the inner 10" as is th e likely 
case for IRAS 12496-7650 (see v an Kempen et al.l(l2006h ) or as- 
sociated with an energetic outflow. 
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Fig. 11. Spectra of '^CO 6-5 over an area of 40" x 40" of Ced 1 10 IRS 4 (left) and BHR 71 ((right). For Ced 1 10 IRS 4, the 
horizontal axes are -10 to 25 km s"' with vertical axes from -4 to 10 K. For BHR 71, horizontal axes are -20 to 10 km s"' with 
vertical axes -5 to 15 K. For BHR 71, the main outflow axis of the red and blue are shown with a dashed and solid lines. 
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Fig. 12. Spectra of '^CO 6-5 over an area of 40" x 40" of IRAS 12496-7650 (left) and RCrA IRS lA((right). RCrA IRS 7A was 
observed using the hexa mode, so the covered area is slightly smaller than that of IRAS 12496-7650 taken in OTF mode. For all 
spectra, the horizontal axes are -15 to 15 km s"' with vertical axes from -5 to 15 K for IRAS 12496-7650 and from -10 to 25 km 
s ' and -5 to 50 K for RCrA IRS 7. For RCrA IRS 7, the main outflow axis of the red and blue lobes are shown with a dashed and 
solid lines. 



5. Outflows 

5.1. Shocks 

iBachiUer & T^faUal ([T999h propose an evolutionary sequence of 
outflows around low-mass protostars, with young deeply embed- 
ded YSOs (Class 0) producing highly collimated and energetic 
outflows, and with more evolved embedded YSOs (Class I) pro- 
ducing outflow s which show les s energetic shocks and a wider 
opening angle dArce et al. I l2007l see also). Observations of the 
shocked '^CO 6-5 and 7-6 gas in the outflow directions confirm 



the scenario that the shocks within the vicinity of the protostar 
grow weaker in energy over time. The Class sources (NGC 
1333 IRAS 2, BHR 71 and RCrA IRS 7) afl show shocked wai-m 
gas in their CO 6-5 and 7-6 lines, in both blue- and red-shifted 
outflow. Of the Class I sources, only HH 46 shows shocked 
gas in its red-shifted outflow. L 1551 IRS 5, one of the most- 
studied molecular o utflows (Moriartv-Schieven & Snell 19881 
IBachiUer etalJll994 . shows little emission in the high-/ lines 
associated with outflow shocks. Although several spectra do 
show a small outflow wing, the integrated emission in the wings 
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Fig. 13. '^CO 6-5 integrated intensity maps of BHR 71 {upper left), Ced 1 10 IRS 4 (upper middle), HH 46 (wpper ng/z/), IRAS 
12496-7650 (/ower Ze/f) and TMR 1 (lower middle). 

Ced t ID IBS 4 fC IJ 2-1 TKR 1 |C I| 3-1 




SO -sa so -sa 

A RA ("] A E!A ("] 



Fig. 14. [C I] 2-1 integrated intensity maps of Ced 1 10 IRS 4 (left) and TMR 1 (ng/if)- Note that the TMR 1 data suffer from a 
de-rotation problem in the calibrator data and are subject to change. 



Table 6. '^CO model predictions 



Source 


J TusdV 


T 

pe 


ik 




Mod. 


Obs. 


Mod. 


Obs. 


TMR 1 


3.6 


4.3 


1.9 


1.9 


IRAS 12496-7650 


7.2 


3.3 


3.6 


1.6 


HH 46 


1.9 


3.9 


1.2 


2.9 


BHR 71 


4.7 


1.4 


2.7 


1.2 


Ced 110 IRS 4 


2.3 


5.2" 


1.2 


1.9 



" Possible outflow contribution, see Fig. [S] and [T3] estimated to be at 
25% of total integrated intensity. 



is not higher than a few cr. All other Class I flows show no sign 
of shocked warm gas. HH 46 and L 1551 IRS 5, even though 
classified as Class I, have massive envelopes of a few Mq, more 



characteristic of Class (|j0rgensen et alj|2002l) . L 1551 IRS 5 is 
believe d to be older and to consist of several s ucessive ejection 
events jBachiller et aI.lll994IWhite et alj|2000h . 

Figure [IT] shows the maximum outflow velocities of the 
'^CO 6-5 and 3-2 lines vs. the bolometric temperature, enve- 
lope mass and outflow force. A clear absence of warm high- 
velocity material is seen for sources with a higher bolometric 
temperature. Only the cold outflow of IRAS 12496-7650 is seen 
at higher Tboi (van Kempen et al. 2006). Similarly, there is also 
a clear relation between the mean outflow force of both red and 
blue outflow lobes, and the maximum velocity seen in both I^CO 
3-2 and '^CO 6-5 emission. Outflow forces are derived from the 
spatial scales and velocities of low-excitation CO line emission 
(1-0 , 2-1 and 3-2) (Ca brit & Bertout 1992, Bourke et al. 199^ 
Ho gerheiide et alJll998i Paper I, van Kempen et al. 2009 sub- 
mitted). 
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Fig. 15. The SEDs and radial profiles at 850 j-im of the sample. Overplotted is the best-fitting envelope model. For Ced 1 10 IRS 4 
and BHR 71, the model spatial 850 fim profile of 1.5 is plotted since data are lacking. Note that this power law index of 1.5 does 
not equal the power law index p. 



5.2. Temperatures of the swept-up gas 

The excitation temperature of the outflowing swept-up gas can 
be derived from the ratios of different CO line wings. As an 
example. Fig. [18] shows the CO 3-2 data from van Kempen et 
al. 2009 submitted and Paper I and CO 6-5 data from this pa- 
per overplotted on the same scales for a few sources. The CO 
6-5 data have not been binned to the larger CO 3-2 beam, so 
the comparison assumes similar volume filling factors of the 
shocked gas; for HH 46 in paper 1 it has been checked that rebin- 
ning gives similar results within the uncertainties. If the density 
is known this excitation temperature can be related to the kinetic 
temperature using the diagnostic plots produce d by the RADEX 
radiative transfer code (I van der Tak et al.ll2007l) . 

As can be seen from the temperature analysis of HH 46 
(Paper I), there are many uncertainties, leading to error bars as 
large as 50 K on the inferred temperatures. The ratios depend on 
the velocity, the optical depth of the line wings and the ambient 
density at different distances from the source. For HH 46, a drop 
in temperature was observed if the density remains constant, but 



it is more plausible that the density is lower at larger radii which 
will result in a constant kinetic temperature with distance. 

Table|2]gives the median temperatures using the extreme ve- 
locities of the line wings with intensities > 3cr, assuming the 
line wings are optically thin and with an am bient cloud density 
of lO"* cm"^. CO 3-2 spectra can be found in Hogerheiide et alj 
( 1998) (TMR 1, L 1551 IR S 5), Knee & Sanden'(2000) (NGC 
1333 I RAS 2'). lPariseet aD(i2006l) (BHR TD. Ivan Kempen etal.l 
(12006') (IRAS 12496-7650), van Kempen et al. 2009 submitted 
(RCrA IRS 7, Ced 1 10 IRS 4) and Paper I (HH 46). 

The main error on these temperature estimates is the optical 
depth in the outflow wings. Deep '^ CO 3-2 observations show 
that outflows can be optically thick. iHogerheiide et alj (Il998h 
find optical depths in the '^CO 3-2 line wings, Twing, on the order 
of 10, while in Paper I T„ing=l-1.8 is found for HH 46. Even if 
the optical depths are assumed to be similar in CO 3-2 and CO 
6-5, kinetic temperatures can be almost 50% lower than those 
given in Table |2l 
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Fig. 17. y„,x of I'CO 6-5 and 3-2 vs. Tboi (left), envelope mass (middle) and outflow force. Triangles represent the blue outflows 
and diamonds the red outflows. Open symbols a re the CO 6-5 lines and filled symbols the CO 3-2 fines. Outflow force from 
ICabrit & Bertou3(ll992h . lHogerfieiide et al.1 (Il998l) (NGC 1333 IRAS 2, L 1551 IRS 5 and TMR D.lBourk e et al.'('1997') (BHR 71) 
and van Kempen e t al. su bmitted (C ed 110IRS4, IRAS 12496-7650, RCrA IRS 7) and lvanKempen et al . (2009) (HH 46). CO 3-2 



from Parise et al. (2006) (BHR 71),'HogerheijdeetalJ ([l998) (TMR 1, L 1551 IRS 5). lKnee & Sandelll (2000) (NGC 1333 IRAS 
2). rvan Kempen et al., (.2006 ) (IRAS 12496-7650) and van Kempen et al. 2009, submitted (HH 46, RCrA IRS 7 and Ced 1 10 IRS 4). 
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Fig. 18. Line wings of CO 6-5 (solid) and CO 3-2 (dashed) for HH 46 (red wing at central position and in the red outflow), Ced 
1 10 IRS 4 (blue at the central position), RCrA IRS 7A (red and blue at the central position). 



Table 7. Outflow kinetic temperature estimates". 



Source 


Ratio* 
K 


rid„(Blue) 
K 


Ratio* 


rk,„(Red) 


RCrA IRS 7 


1 


>200 


2.5 


160 


NGC 1333 IRAS 2 










L 1551 IRS 5' 


>8 


<50 


9 


40 


TMR 1'' 


1.5 


180 


>4 


120 


IRAS 12496-7650 


>8 


<65 






HH 46'' 


6 


100 


3 


140 


BHR 71 










Ced 110 IRS 4 


3 


140 


>7 


<80 



" Estimated with an ambient density of 10"* cm"''. 
* CO 3-2/6-5 line wing ratio . Ratios are averaged ove r velocity and 
distance from the source. See Ivan Kempen et alT ( l2009l) for more dis- 
cussion. 

' CO 4-3/6-5 ratio. 

'' Opt ical depth of 8 inferred for '-CO 3-2 wings jHogerheiide et al] 

Il998h 

" See lvan Kempen et alj ( 1200^ . 



For a more thorough discussion about the density of the sur- 
rounding envelope and cloud material, see Paper I. Our choice of 



the density of 10 cm is based on th e typic al densities found 
in model envelopes of |j0rgensen et al.l (l2002h at distances of a 
few thousand AU, corresponding to the 10 K radius where the 
envelope merges with the surrounding molecular cloud. Higher 
densities will produce lower outflow temperatures (see paper I 
for more extensive discussion. 



A theoretical study by iHatchell et al.l ( 1199 9') using a jet- 
driven bow shock model suggests that the kinetic temperatures 
of molecular outflows are typically of order of 50-150 K along 
the axis and rising toward the location of the bow shock, with 
overall values increasing with higher jet velocities. The temper- 
atures found in Table [T] agree well with these predicted temper- 
atures. For several outflows, e.g., IRAS 12496-7650 (blue) or 
Ced 110 IRS 4 (red), temperatures are significantly lower, how- 
ever. The lowest outflow temperatures are found for the flows of 
L 1551 IRS 5, where CO 4-3/6-5 ratios are on the order of 8 or 
higher. At an assumed density of 10"* cm"-', this corresponds to 
kinetic temperatures of 50 K and lower. If both line wings are 
optically thin, densities must be lower than 10^ cm"^ to produce 
temperatures of ~100 K, observed for other flows. 
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6. Heating processes in the molecular outflow and 
protostellar envelope 

As discussed in § 4, several sources show quiescent, narrow 
'^CO 6-5 and 7-6 emission that is more intense than can be 
produced by an envelope model. Moreover, strong narrow high- 
J CO emission is observed ofF-source along outflow axes for 
most sources. In Paper I, we proposed that for HH 46 the qui- 
escent narrow '^CO 6-5 and 7-6 line emission originates within 
the outflow cavity walls heated to 250^00 K by an enhance- 
ment factor Go with respect to the standard interstellar radia- 
tion field of a few hundred. The UV photons can be created by 
jet shocks in the outflow cavities as well as in the disk-star ac- 
cretion boundary layer near the central protostar . This heating 
method was first proposed bv lSpaans et al.l ( 1199 5'). but extended 
by Paper 1 to both the inner envelope as well as the outflow cav- 
ity walls much further from the central star The data presented 
in this paper show that photon heating is present in other proto- 
stars as well, especially in outflow cavities. Even in sources with 
little to no outflow, such as Ced 110 IRS 4, relatively strong nar- 
row '^CO 6-5 and 7-6 lines are seen at positions not associated 
with the central protostar, see Fig. [161 The presence of signifi- 
cant radiative 'feedback' from the protostar on its surroundings 
may have consequences for the collapse of the envelope, limit- 
ing the accretion ra te and mass o f the star and suppressing disk 
fragmentation (e.g.. lOffner et al.ll2009 ). 

The origin of such quiescent high-/ CO emission at the 
source velocity is physically different from both the ther- 
mal emission of the protostellar envelope (see § 6.4), emis- 
sion from shocks present at the w orking surfaces of outflows 
jReipurth & Ragalll999llRaga et alj|2007l) . 

Slow (vs - 5 - 10 km s~^) C-shocks may produce sim- 
ilar q uiescent levels of CO 6-5 and 7-6 emission in out- 
flows jDraine & Robergd [T9841 ISpaans et al.lil995l) . However, 
the presence of quiescent ^^CO 6-5 and 7-6 emission in the 
envelopes of TMR 1 and Ced 110 IRS 4, both of which show 
little to no spatially resolved outflow emission in the CO 3-2 
line (Hogerheiide et al. 1998, and van Kempen et al. submit- 
ted), is more easily explained with the photon heating scenario 
than with slow C-shocks. In addition, the narrow line widths for 
other sources argue against this scenario. 



6. 1. Envelope and outflow of BHR 71 



Fig. 16. The contribution of the envelope to CO lines at the 
source position with 7up ranging from 1 to 19. Integrated emis- 
sion is shown with dashed and solid lines. Dashed lines are mod- 
elled lines which show excessive self-absorption. Solid lines rep- 
resent gaussian fitted to the Une wings of the envelope model 
and thus represent the strict upper limit to the CO line emission 
from the protostellar envelope. Beam-sizes were taken to be 10" 
for all transitions except for J^p <3, for which a beam of 20" 
was used. Stars indicate the integrated quiescent emission of the 
lines observed with CHAMP^ from Table [3] at the source posi- 
tion, corrected for red- and blue-shifted emission, while triangles 
are hterature data (see Table |2|. 



The proposed model that photon heating takes place both in the 
outflow cavities and the inner envelope is clearly illustrated by 
the observations of BHR 71. At the north and south position 
of the outflow, shocked gas is clearly detected, but almost no 
shocked gas is seen near the central star Quiescent gas with 
Ay=1.5 km s ' is observed both within the inner envelope and 
along the outflow. In the '^CO 7-6 map, the photon heating in 
the envelope can be strongly identified by the central contour 
BHR 71 is the only source for which the quiescent emission 
at the outflow position is (much) brighter than that seen at the 
position of the envelope. It is likely that strong shocks produc- 
ing copious UV are present in the main outflow. This outflow 
dBourke et al.ll 19971) is clearly detected at the edges of the map, 
but the secondar y weaker outflow, a ssociated with the 1RS2 po- 
sition, as seen bv' Parise et al.l (12006') is not detected down to the 
noise levels. Shocked CO 6-5 and 7-6 emission is also only seen 
at relatively low levels in the inner 20" and strongly increases at 
40" away from the source for both outflows. 
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6.2. The 'fossil' outflow of L 1551 IRS 5 

The outflow of L 1551 IRS 5 has been considered 
an example of an older outfl o w, due to its larg e size 
dMoriartv-Sch ieven & SnelT '1988', 'Fridl und et all Il989l) . large 
opening angle (Bachiller et al. 1994) an d other subm i llime- 
ter properties C abrit & Bertout 1992, Bach iUer et al] 11994 
iHogerheiide et aTl i 1998). White et al, (2000) constructed a de- 
tailed model from a wide variety of space- and ground-based 
IR continuum and spectroscopic observations, including a wide 
opening angle of the outflow (~50°) and low densities inside the 
cavity. This view is clearly confirmed by the low derived temper- 
atures of 40 to 50 K of the shocked gas, which is much cooler 
than in other sources. Little high-7 CO emission is seen at the 
off-positions, and the maps are dominated by the emission of 
the central source. Comparison of the peak temperatures of CO 
4-3 dHogerheiide et alJll998h . 6-5 and 7-6 shows that Tpeak at 
high-7 transitions is similar to that at low-/ CO transitions. This 
suggest that even though the outflow is present, the density in 
the surrounding cloud must be quite low. 



6.3. The PDR of RCrA IRS 7 

The lines of RCrA IRS 7 are an order of magnitude stronger 
than those of equivalent sources in different clouds. Even at off"- 
positions at the edge of the observed area, integrated intensities 
larger than 300 K km s ' are seen. Even with the high luminosity 
of ~20-30 L0, these integrated intensities cannot originate from 
heating by RCrA IRS 7 itself. A likely explanation is the prox- 
imity of the A5 star RCrA (Lboi = 130 Lq) at 36 " (4,500 AU). 
RCrA was observed by ISO-LWS, and many strong high-7 CO 
lines (CO 14-13 to 2 1-20) wer e detected in the large beam of 
ISO dLorenzetti et al.i ri999. Gia nnini et al., 1999). Models were 
put forward with the emission originating in relatively small 
(-0.002 pc), dense (>10'' cm -^) and hot (T > 300 - 1000 K) 
regions (iGiannini et al.|[T999l) . However, the spatial distribution 
of the CO 6-5 and 7-6 seen in Fig.[8]does not agree with this hy- 
pothesis. A much more likely explanation is that RCrA itself ir- 
radiates the outer edges of the cloud and envelope around RCrA 
IRS 7a, creating a PDR at its surf ace, much like th e case of in- 
termediate mass sources in Orion (|j0rgensen et alj|2006,) . 



6.4. Presence of [C I] 2-1 

A limit on the amount of FUV/X-ray emission that is available 
to dissociate CO and H2 can be derived from the presence of the 
[C I] 2-1. Photodissociation of CO can only occur at 912-1100 
A dvan Disho eck & Black 1988), so the absence of strong [C I] 
2-1 emission in the outflow cavities of most protostars suggests 
that the radiation field in outflows does not produce sufficiently 
energetic radiation, but still heats the cavity walls to a few hun- 
dred K. This would also limit the shock velocities to <90 km 
s \ since hi gher shock velocities produ ce CO and H2 dissocat- 
ing photons dNeufeld & Dalgarnolll989l) . 

[C I] is detected in the inner 10" for HH 46, TMR 1 and 
Ced 110 IRS 4 at levels of 2 K km s '. For IRAS 12496-7650 
and BHR 71, no [CI] 2-1 was detected. The low emission in all 
sources can be accounted for by a C abundance of a few times 
10"^ with respect to H2. Such abundances can be m aintainted by 
low U V levels produced by cosmic ray radiation dFlower et alj 
fT994l) . 



7. Conclusions 

In this paper, we have presented the first ^^CO 6-5 and 7-6 maps 
of a sample of 8 low-mass protostars with a large range of lumi- 
nosities, evolution and densities, as well as several isotopologue 
and [C I] observations. All observations have been carried out 
with the CHAMP+ instrument. The main conclusions of this pa- 
per are 

- Warm gas, as traced by '^CO 6-5 and 7-6, is present in 
all observed protostars at the central position. Three diff'er- 
ent origins of the warm gas emission are found: (i) the in- 
ner envelope heated passively by the protostellar luminos- 
ity; (ii) shocked gas in the outflow; and (iii) quiescent gas 
heated by UV photons. This latter mechanism, as detailed in 
ISpaans et al.l dl995l) and in Paper I, generally dominates the 
extended high-7 CO emission. 

- Envelope models show that for HH 46 and Ced 1 10 IRS 4, 
passive heating of the envelope is insufficient to explain the 
observed '-CO 6-5, 7-6 and '^CO 6-5 lines, requiring heat- 
ing of the envelope by UV photons even at the (0,0) position. 

- Photon heating of the cavity walls takes place on arcmin 
scales for the outflow cavities of several other Class and 
Class I protostars, as seen at positions off-source of BHR 71, 
NGC 1333 IRS 2, HH 46 and L 1551 IRS 5. The necessai-y 
UV photons are created by internal jet shocks and the bow 
shock where the jet interacts with the ambient medium, in 
addition to the disk-star accretion boundary layer. The distri- 
bution of the quiescent CO 6-5 and 7-6 emission of BHR 7 1 , 
with narrow emission stronger at larger distances from the 
source, confirms the hypothesis that UV photons necessary 
for the heating can originate both mechanisms. This suggests 
that photon heating is present in all outflows. 

- The lack of [C I] 2-1 emission in the outflows constrains 
the production of energetic CO dissociative photons in the 
shocks. The observed [C I] emission at the source position 
can be accounted for by a low atomic C abundance that is 
maintained by cosmic ray induced UV photodissocation of 
CO. 

- Shocked '^CO 6-5 and 7-6 gas only exists in large quan- 
tities in flows of more massive sources (NGC 1333 IRAS 
2, BHR 71, HH 46 and RCrA IRS 7), while low-mass 
flows do not show much shocked gas that emits in the 
'•^CO 6-5 and 7-6 lines. Within the proposed model of 
iBachiller & TafaUal ([T999) . the outflows of Class I are more 
evolved and are driven with much less energy, producing 
much weaker shocks. This is also reflected in the decreasing 
maximum velocity of '^CO 6-5 with lower outflow force. 

- Kinetic temperatures of ~100 K are found for the molecu- 
lar gas in most flows studied here. Such temperatures agree 
with expected condi tions of jet driven flows modelled by 
^Hat chell et ai][T999l) . Only the flows of L 1551 IRS 5 and 
IRAS 12496-7650 are much colder (<50 K), in agreement 
with the 'fossil', empty nature of the flows. 

- The very strong intensities at all positions of RCrA IRS 7, an 
order of magnitude higher than can be produced by passive 
heating, must be due to a significant PDR region near the 
source. Likely, RCrA itself heats the outside of the RCrA 
IRS 7 envelope and cloud region. 

Spatially resolved observations with APEX-CHAMP+ of 
high-7 CO and isotopologues clearly provide unique insight 
into the structure and physical processes of low-mass proto- 
stars, as well as their feedback on the surroundings. The ubiq- 
uitous photon heating of the gas found here may have some 
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consequences for the amount of envelope gas that can eventu- 
ally collapse. Detailed multi-dimensional radiative transfer mod- 
elling, including outflow cavities are needed to directly adress 
this question . Such models are n ow being developed for high- 
mass YSOs (IBruderer et al.ll2009l) . and will be extended to low- 
mass YSOs. The UV radiation will also significantly affect the 
chemistry along the outflow walls, dissociating molecules like 
H2O to OH and O and HCN to CN. In the future, the Herschel 
Space Observatory will be able to map even higher J CO lines 
as well as H2O, OH and [O 1] lines in the far-lR between 60 and 
600 /zm with the PACS and HIFI instruments. Together, radia- 
tive transfer modelling, these data sets and future interferomet- 
ric observations with ALMA will be able to fully characterize 
the physical and chemical structure of the earliest, deeply em- 
bedded stages of star formation. 

Acknowledgements. TvK and astrochemistry at Leiden Observatory are sup- 
ported by a Spinoza prize and by NWO grant 614.041.004. The staff at APEX 
and many members of the Bonn submillimeter group are thanked for the ex- 
cellent support during observations. CHAMP^ is constructed with funds from 
NWO grant 600.063.310.10. We thank Ronald Stark for continued support of 
the CHAMP+ project. 



References 

Adams, F. C, Lada, C. J., & Shu, F. H. 1987, ApJ, 312, 788 
Arce, H. G. & Sargent, A. I. 2005, ApJ, 624, 232 

Arce, H. G., Shepherd, D.. Gueth, F., et al. 2007, in Protostars and Planets V, ed. 

B. Reipurth, D. Jewitt, & K. Keil, 245-260 
Bachiller, R. 1996, ARA&A, 34, 111 

Bachiller, R. & Tafalla, M. 1999, in The Origin of Stars and Planetary Systems, 

ed. C. J. Lada, , & N. D. Kylafis (Kluwer, Dordrecht), 227 
Bachiller, R., Tafalla, M., & Cemicharo, J. 1994, ApJ, 425, L93 
Bisschop, S. E., J0rgensen, J. K., van Dishoeck, E. F, Sl de Wachter, E. B. M. 

2007, A&A, 465,913 

Blake, G. A., Sandell, G., van Dishoeck, E. F, et al. 1995, ApJ, 441, 689 
Boogert, A. C. A., Hogerheijde, M. R., Ceccarelli, C, et al. 2002, ApJ, 570. 708 
Bottinelli, S., Ceccarelli, C, Lefloch, B., et al. 2004, ApJ, 615, 354 
Bottinelli, S., Ceccarelli, C, Wilhams, J. R, & Lefloch, B. 2007, A&A, 463, 601 
Bourke, T. L. 2001, ApJ, 554, L91 

Bourke, T. L., Gai'ay, G., Lehtinen, K. K., et al. 1997, ApJ, 476, 781 
Brinch, C, Crapsi, A., J0rgensen, J. K., Hogerheijde, M. R., & Hill, T. 2007, 
A&A, 475,915 

Bruderer, S., Benz, A. O., Doty, S. D., van Dishoeck, E. F. & Bourke. T. L. 

2009, ApJ, 700, 872 
Butner, H. M., Evans, II, N. J., Lester, D. R, Levreault, R. M.. & Strom, S. E. 

1991, ApJ, 376, 636 
Cabrit, S. & Bertout, C. 1992, A&A, 261, 274 

Ceccarelli, C, Boogert, A. C. A., Tielens, A. G. G. M., et al. 2002, A&A, 395, 
863 

CeccarelU, C, Castets, A., Caux, E., et al. 2000, A&A, 355, 1129 
Ceccarelh, C, Caux, E., White, G. J., et al. 1998, A&A, 331, 372 
Chandler, C. J. & Richer, J. S. 2000, ApJ, 530, 851 

Crapsi, A., van Dishoeck, E. R, Hogerheijde, M. R., Pontoppidan, K. M., & 

Dullemond, C. P 2008, Ai'Xiv e-prints, 801 
Di Francesco, J., Johnstone, D., Kirk, H., MacKenzie. T, & Ledwosinska, E. 

2008, ApJS, 175, 277 

Draine, B. T. & Roberge, W. G. 1984, ApJ, 282, 491 

Dunham, M. M., Crapsi, A.. Evans, II, N. J., et al. 2008, ArXiv e-prints, 806 
Enoch, M. L., Young, K. E., Glenn, J., et al. 2006, ApJ, 638, 293 
Evans, N. J., Harvey, P. M., Dunham, M. M., et al. 2007, c2d delivery document 
Flower, D. R., Le Bourlot, J., Pineau Des Forets, G., & Rouefl", E. 1994, A&A, 
282, 225 

Fridlund, C. V. M., Sandqvist, A., Nordh, H. L., & Olofsson, G. 1989, A&A, 
213, 310 

Froebrich, D. 2005, ApJS, 156, 169 

Furlan, E., McClure, M., Calvet, N., et al. 2008, ApJS, 176, 184 
Giannini, T, Lorenzetti, D., Tommasi, E., et al. 1999, A&A, 346, 617 
Giannini, T, Nisini, B., & Lorenzetti, D. 2001, ApJ, 555, 40 
Groppi, C. E., Hunter, T. R., Blundell, R., & Sandell, G. 2007, ApJ, 670, 489 
Glisten, R., Baryshev, A., Bell, A., et al. 2008, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 7020, Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference Series 
Gutermuth, R. A., Myers, P C, Megeath, S. T, et al. 2008, ApJ, 674, 336 



Hatchell, J., Fuller, G. A., & Ladd, E. F 1999, A&A, 344, 687 
Hatchell, J., Fuller, G. A., Richer, J. S., Harries, T. J., & Ladd, E. F 2007, A&A, 
468, 1009 

Heathcote, S., Morse, J. A., Hartigan, P, et al. 1996, AJ, 112, 1141 
Henning, T. & Launhardt, R. 1998, A&A, 338, 223 
Henning, T, Pfau, W., Zinnecker, H., & Prusti, T. 1993, A&A, 276, 129 
Hiramatsu, M., Hayakawa, T, Tatematsu, K., et al. 2007, ApJ, 664, 964 
Hirano, N., Liu, S.-Y., Shang, H., et al. 2006, ApL 636, L141 
Hogerheijde, M. R. & van der Tak, F F. S. 2000, A&A, 362, 697 
Hogerheijde, M. R., van Dishoeck, E. F, Blake, G. A., & van Langevelde, H. J. 
1998, ApJ, 502,315 

Hogerheijde, M. R., van Dishoeck, E. F, Salverda. J. M., & Blake, G. A. 1999, 

ApJ, 513, 350 
Ivezic, Z. & Elitzur, M. 1997, MNRAS, 287, 799 
J0rgensen, J. K. 2004, A&A, 424, 589 

J0rgensen, J. K., Johnstone. D., van Dishoeck, E. F, & Doty, S. D. 2006, A&A, 
449, 609 

J0rgensen, J. K., Schoier, F L., & van Dishoeck, E. F 2002, A&A, 389, 908 
J0rgensen, J. K., Schoier, F L., & van Dishoeck, E. F 2004, A&A, 416, 603 
J0rgensen, J. K., Schoier, F L., & van Dishoeck, E. F 2005, A&A, 435, 177 
Kasemann, C, Glisten, R., Heyminck, S., et al. 2006, Presented at the Society of 

Photo-Optical Instrumentation Engineers (SPIE) Conference, 6275 
Kenyon, S. J., Gomez, M., Marzke, R. O., & Hartmann, L. 1994, AJ, 108, 251 
Knee, L. B. G. & SandeU, G. 2000, A&A, 361, 671 
Lee, C.-R & Ho, R T. P 2005, ApJ, 624, 841 

Lee, J.-E., Di Francesco, J., Bourke, T. L., Evans, II, N. J., & Wu, J. 2007, ApJ, 
671, 1748 

Lehtinen, K., Haikala, L. K., Mattila, K., & Lemke, D. 2001, A&A, 367, 311 
Lehtinen, K., Harju, J., Kontinen, S., & Higdon, J. L. 2003, A&A, 401, 1017 
Liu, M. C, Graham, J. R., Ghez, A. M., et al. 1996, ApJ, 461, 334 
Lorenzetti, D., Tommasi, E., Giannini, T, et al. 1999, A&A, 346, 604 
Luhman, K. L., Allen, L. E., Allen, P R., et al. 2008, ApJ, 675, 1375 
Maret, S., Ceccarelh, C, Caux, E., et al. 2004, A&A, 416, 577 
Moriarty-Schieven, G. H. & Snell, R. L. 1988, ApJ, 332, 364 
Motte, R & Andre, P 2001, A&A, 365, 440 
Neufeld, D. A. & Dalgarno, A. 1989, ApJ, 340, 869 
Nisini, B., Benedettini, M., Giannini, T, et al. 1999, A&A, 350, 529 
Nisini, B., Giannini, T, & Lorenzetti, D. 2002. ApJ, 574, 246 
Nutter, D. J., Ward-Thompson, D., & Andre, R 2005, MNRAS, 357, 975 
Offner. S. S. R., Klein, R. I., McKee, C. R, & Krumholz. M. R. 2009, ArXiv 
e-prints 

Osorio, M., D'Alessio, P., MuzeroUe, J., Calvet, N., & Hartmann, L. 2003, ApJ, 
586, 1148 

Ossenkopf V. & Henning, T. 1994. A&A, 291, 943 
Parise, B., Belloche, A., Leurini, S., et al. 2006, A&A, 454, L79 
Raga, A. C, de Colle, R, Kajdic, R, Esquivel, A., & Canto, J. 2007, A&A, 465, 
879 

Reipurth, B., Chini, R., Krugel, E., Kreysa, E., & Sievers, A. 1993, A&A, 273, 
221 

Reipurth, B. & Raga. A. C. 1999, in NATO ASIC Proc. 540: The Origin of Stars 
and Planetary Systems, ed. C. J. Lada & N. D. Kylafis, 267-+ 

Reipurth, B., Rodriguez, L. R, Anglada, G., & Bally, J. 2002, AJ, 124, 1045 

Robitaille, T. P., Whitney, B. A., Indebetouw. R.. Wood, K.. & Denzmore, P. 
2006, ApJS, 167, 256 

Sandell, G. & Knee, L. B. G. 2001, ApJ, 546, L49 

Schoier, R. L., J0rgensen, J. K., Pontoppidan. K. M., & Lundgren, A. A. 2006, 
A&A, 454, L67 

Schoier, R L., J0rgensen, J. K., van Dishoeck, E. R, & Blake, G. A. 2002, A&A, 
390, 1001 

Schoier, R L.. van der Tak, R. R S., van Dishoeck, E. R, & Black, J. H. 2005, 
A&A, 432, 369 

Schuster, K. R, Harris, A. I.. Anderson, N., & Russefl, A. P G. 1993, ApJ, 412, 
L67 

Schuster, K.-R, Russefl, A. P G., & Harris, A. I. 1995, Ap&SS, 224, 117 
Shirley, Y. L., Evans, N. J., Rawlings, J. M. C, & Gregersen, E. M. 2000, ApJS, 
131, 249 

Shirley, Y. L., Evans, II, N. J., & Rawhngs, J. M. C. 2002, ApJ, 575, 337 
Spaans, M., Hogerheijde, M. R., Mundy, L. G., & van Dishoeck, E. F. 1995, 

ApJ, 455, L167+ 
Stark, R., Sandell, G., Beck, S. C, et al. 2004, ApJ, 608, 341 
Stauber, P., Doty, S. D., van Dishoeck, E. F, J0rgensen, J. K., & Benz, A. O. 

2004, A&A, 425, 577 
Terebey, S., Chandler, C. J., & Andre, R 1993, ApJ, 414, 759 
van der Tak, F. F. S., Black, J. H., Schoier, F. L., Jansen, D. J., & van Dishoeck, 

E. F 2007, A&A, 468, 627 
van Dishoeck, E. F & Black, J. H. 1988, ApJ, 334, 771 

van Dishoeck, E. F, Blake, G. A., Jansen, D. J., & Groesbeck, T. D. 1995, ApJ, 
447, 760 



T.A. van Kempen et al.: The warm gas within young low-mass protostars 



van Kempen, T. A., Hogerheijde, M. R., van Dishoeck, E. F., et al. 2006, A&A, 

454, L75 

van Kempen, T. A., van Dishoeck, E. F., Salter, D. M., et al. 2009, A&A, 498, 

167 

White, G. J., Liseau, R., Men'shchikov, A. B., et al. 2000, A&A, 364, 741 
Whitney, B. A., Wood, K., Bjorkman, J. E., & Cohen, M. 2003a, ApJ, 598, 1079 
Whitney, B. A., Wood, K., Bjorkman, J. E., & Wolff, M. J. 2003b, ApJ, 591, 
1049 




VelpciLy (kasi/a) 



